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Abstract 

In high energy hadronic colhsions a scalar or pseudoscalar Higgs boson, ^ = H, 
A, can be efficiently produced via gluon fusion, which is mediated by heavy quark 
loops. We here consider double real emission corrections to <I> = ^ production, 
which lead to a Higgs plus two-jet final state, at order . Full quark mass effects 
are considered in the calculation of scattering amplitudes for the CP-odd Higgs boson 
A, as induced by quark triangle-, box- and pentagon-diagrams. They complement 
the analogous results for a CP-even Higgs boson H in Ref. [T]. Interference effects 
between loops with top and bottom quarks as well as between CP-even and CP-odd 
couplings of the heavy quarks are fully taken into account. 
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1 Introduction 



One of the prime tasks of the CERN Large Hadron Colhder (LHC) is the search for the 
origins of the spontaneous breaking of the electroweak SU (2) x U{1) gauge symmetry and, 
once such particles are found, the study of one or several Higgs bosons as the remnants of 
the symmetry breaking mechanism. Among the various Higgs boson production channels 
the gluon fusion and the weak boson fusion processes have emerged as the most promising 
channels for Higgs boson discovery at the LHC [21 El H], and they are equally valuable 
for the study of its properties, like the measurement of its couplings to gauge boson and 
fermions [5l [6] . 

In weak boson fusion (WBF), qq^qqH mediated by t-channel W or Z exchange, the 
two forward tagging jets arising from the scattered quarks provide a tell-tale signature 
which can be used for efficient background rejection [7j. The same Hjj signature can also 
arise in gluon fusion events, via O(a^) real emission corrections to gg^H which, within 
the Standard Model (SM), is mediated mainly by a top quark loop. For a Higgs boson 
which is lighter than the top quark, the resulting Hjj cross section can be determined to 
good approximation by an effective Lagrangian of energy dimension five, which is given 
by Pi 

where G^j, denotes the gluon field strength and G'*^^ = 1/2 0^^6'^'^^°' its dual. The two 
terms result from a ttH and a ti'j^tA coupling of the (pseudo) scalar Higgs, respectively, 
and they lead to distinctively different distributions of the azimuthal angle between the 
two jets: the CP-even Hgg couphng produces a minimum for = ±7r/2 while a CP-odd 
Agg coupling leads to minima at (pjj = and ivr. These distinctions become important 
in two-Higgs-doublet models (2HDM) like the minimal supersymmetric extension of the 
standard model (MSSM), where a CP-odd Higgs, A, appears in addition to a light and a 
heavy, neutral CP-even Higgs, h and H: the azimuthal angle distribution of $jj events 
allow to differentiate between a CP-even Higgs ($ = h, H) or a CP-odd one ($ = A). 

For large Higgs boson masses {rriH ^ rrit) the full quark mass dependence of the loop 
diagrams must be calculated for reliable predictions, and the same is true for large ratios of 
the two vacuum expectation values, Vu/vd = tan /3, where bottom quark loops provide the 
dominant contributions to qq^qqH, qg^qgH and gg^ggH amphtudes. For a CP-even 
Higgs boson, these calculations were performed in Ref. [1]. The purpose of the present 
paper is to present the corresponding results for a CP-odd Higgs boson, $ = A, or, more 
precisely for an underlying Higgs coupling to quarks derived from the Yukawa Lagrangian 
C = Uqqi'-fr^qA. By combining the present results with those for a CP-even Higgs, the 
quark loop induced contributions to $jj production can be calculated for an arbitrary 
Yukawa coupling of the form 

Yukawa = Q {Vq + hbVg) (1-2) 

Our results are implemented in a parton level Monte Carlo program which is part of 
the VBFNLO program package [10]. This numerical implementation allows to calculate 
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$jj production cross sections in hadronic collisions including top- and bottom-quark loop 
contributions for arbitrary combinations of the Yukawa couplings Uq and ijq {q = t,b). 

Our paper is organized as follows. In Section [2] we first define the models in which 
we consider pseudoscalar Higgs production. We then provide an outline of the calcu- 
lation of the scattering amplitudes for the three basic subprocesses, qq^qqA, qg^qgA 
and gg^ggA. Further details on the various loop contributions are relegated to the Ap- 
pendices. We have performed a number of analytic and numerical consistency checks 
on our calculation: they are described in Section [31 The main phenomenological results 
are presented in Section HI for pp scattering at the LHC with a center of mass energy 
of -^i = 14 TeV. For various combinations of top- and bottom quark contributions, 
parameterized by tan/3, we provide integrated Ajj cross sections but also differential dis- 
tributions. Results are also presented for general ^jj events, i.e. for the production of a 
Higgs boson with arbitrary CP-violating couplings to the third generation quarks. Final 
conclusions are drawn in Section 



2 Outline of the calculation and matrix elements 

The production of the CP-odd Higgs boson A in association with two jets, at order 
aj, proceeds in analogy to the production of the CP-even Higgs boson Hsm of the SM. 
The HsmJJ production processes with full heavy quark mass effects were considered in 
Ref. [1], and we here closely follow the framework and the notation introduced there. We 
consider the production subprocesses 

qq -> qqA , qQ qQA , qg qgA , gg ggA , (2.1) 

and all crossing-related processes. Here the first two entries denote scattering of identical 
and non-identical quark flavors. The Higgs boson A is produced by massive quark loops, 
for which only the third quark generation is taken into account. Furthermore, within the 
MSSM, massive squark loops can safely be neglected, because their contribution sums to 
zero at amplitude level in the production of a CP-odd Higgs boson [TT]. In the 2HDM, 
up- and down-type quark Yukawa couplings depend on the ratio of vacuum expectation 
values, tan/3 = Vu/vd, via the relations 

1. 2HDM of type I: 

~i cot /3 A ~i ^ n\ 

yA,uu = yA,dd = — —^d , (2.2) 

2. 2HDM of type II (MSSM): 

~ii cot/3 tan/3 

yA,uu = and = "^d • (2-3) 



In the 2HDM of type I, Yukawa couplings for up-type and down-type quarks are suppressed 
equally at large tan /3 compared to the 2HDM of type II, where only the up-type Yukawa 
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coupling is suppressed but the down-type Yukawa coupling is enhanced. Due to this 
enhancement, loops with bottom quarks can also provide significant contributions to cross 
sections. 

In the calculation of the subprocesses listed in fl2.ip . three different loop topologies 
appear: the triangle-, box-, and pentagon diagrams of Fig. [T] The contributing Feynman 




Figure 1: Examples of Feynman graphs contributing to A + 2 jet production via 
gluon fusion. 



diagrams can be easily built from the simpler QCD dijet processes at leading order. One 
needs to insert the Higgs-gluon triangles into the gluon propagators of the 2^2 tree-level 
diagrams in all possible ways or one replaces a triple gluon or four gluon vertex by box 
or pentagon graphs in all possible ways. Charge-conjugation related diagrams, where the 
loop momentum is running clockwise and counter-clockwise, can be counted as one by 
exploiting Furry's theorem [12]. This effectively reduces the number of diagrams by a 
factor of two. Furthermore all diagrams are UV-convergent and, due to the finite quark 
mass in the loops, also IR-convergent. All coupling constants and loop factors which 
appear can conveniently be absorbed into an overall factor 

Ff = ^mfhf^ = Amfhja',, (2.4) 

where f = b,t labels the heavy quark fiavor of a particular loop. In the following, we use 
the MSSM couplings of Eq. (12. 3p . i.e. we set ht = cot /3 mt/v and hb = tan/3 rrib/v. By 
pulling out a loop factor Amf/lQir"^, we anticipate that the Dirac trace of all loops requires 
a quark mass insertion to compensate the helicity fiip induced by the $// coupling. 

2.1 Subprocesses qQ—^qQA and qq^qqA 

The subprocess qQ^qQA, depicted in Fig. (U^a), is the simplest contribution to A + 2 jet 
production. Following Ref. [1], the amplitude for different flavors can be written as 

= J2 Ffj,'lJ^iT,,,M,Q2,mj)tl,Jl,^ = A^^tl^^tl,^ . (2.5) 

f=t,b 

Using the notation and formalism for the spinor algebra of [I3], the external quark lines 
can be expressed by effective quark currents and 743^ as given in [1]. The triangle 



4 



tensor T'^^^^{qi,q2,mf) (see Appendix Fig. |9]) has the simple form 

T^^^^'iq,,q2,mf) = e'^™'^' dq^.q^^nif) . (2.6) 

Here, Cq denotes the scalar three-point function and e'^i^'2?ig2 jg ^y^q totally anti-symmetric 
tensor (Levi-Civita symbol) in four dimensions contracted with attached gluon momenta 
qi and q2. The = A° /2 are color generators in the fundamental representation of 
SU(A^), = 3 and the overall factor 

Ff = S, S2 Ss S, A^pIp'^pI Ff (2.7) 

includes normalization factors of external quark spinors. Here, the denote physical 
momenta describing phase space and the wave functions of fermions and bosons while Pi 
is used for momenta appearing in the momentum flow in Feynman diagrams. Both sets 
of momenta are related by the sign factors Si [13] 

P^ = S,p,, (2.8) 

with Si = +1 for fermions and Si = —1 for anti- fermions. The factor Ff is given in 
Eq. (12. 4p . For identical quark flavors, one has to keep in mind Pauli interference 

Aqq — All fa .a _ Aqq .a .a fc^ qN 

The squared amplitude, summed over initial- and final-particle color, becomes 

\^''? = (|^2143P + |A123P)^^^ + 2Re(^2143^:i23)^^^ • (2-10) 

color 



2.2 Subprocess qg—^qgA 



Polarization vectors of external gluon lines with a triangle insertion can be expressed by 
effective polarization vectors 

1 



Co{qi, -{qi + P),mf) 



(2.11) 



{q^ + Py 

which replace the polarization vectors ef of the underlying 2— !'2 process for gluons i = 1,2. 
Here qi is the external gluon momentum while P denotes the momentum of the Higgs 
boson. The expression for the amplitude of graphs with a triangle insertion adjacent to a 
three-gluon vertex differs slightly from that in [1] due to the emergence of the Levi-Civita 
symbol 

= Pf [^"'' 1^ [^1^ ■ ^2 J21 ■ q2 - eiA • J21 £2 ■ (p2 - Pi) - eiA ■ q2 J21 ■ £2 



- 2 



+ 2 F 7^^ 



e2A ■ ei J21 ■ qi - e2A ■ qi J21 ■ ei - 62^1 ■ J21 ei ■ {P2 - Pi] 

ei ■ 62 q^'q^'{p2 - PiY' + {t2 ■ qi ef - ei ■ ) 
Cq{p2 -pi,qi + q2,mf) 



X iqi + q2T'{p2-pi] 



M4 



{qi + 92) 



(2.12) 
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Further expressions for amphtudes of graphs with a triangle insertion can be taken from [T] 
replacing e^^ by e^^(mj). The tensor structure of the box diagram in Fig. [H^b) is given 
by B^j^fj_^^.^{qi,q2,q3,mf), which is shown pictorially in Fig. [TU] and given explicitly in 
Appendix [Bl Finally the color structure of the qg-^qgA amplitude is given by 



A''^ = {en'"') .^^^Al^ + {f'H''') .^.^All with A'^' 



/ 



19 
f 



The indices 12 and 21 label amplitudes with interchanged external gluons. 
resulting color-summed squared amplitude takes the form 

,2 



EM' 

color 



■^12 



■^21 



{N^ - ly 



2Re 



A19( A19Y 



4N 



(2.13) 



Thus, the 



(2.14) 



2.3 Subprocess gg—^ggA 

After inserting suitable loop topologies and application of Furry's theorem, this process 
contains 19 graphs with triangle insertions, 18 box contributions and 12 pentagon dia- 
grams. The pentagon diagrams Fig. [D^c) enter via the pt^^i^^tJ.aiJ.A tensor (see Appendix 
Fig. [H]). Full expressions and diagrams can be looked up in [16]. The contributing color 
structures to the process gg-^ggA can be expressed by the real-valued color coefficients 
Ci defined in [T] 

C2 = tT[t^H''H''*e^] + tT[t^H^H''H''^] , (2.15) 
cg = tT[t''H''H''H'''] +tr[ru"3t"2t"^] . 
Evaluation of the color traces yields 

Cl 



C3 



4 




_ 1^0-1 0-2 '^03 14' _ 










_ ^aiasl ^a4a2l 






4 


iV " 


_ j^aiail ^a2a-il 









caia4l i'a20,zl 



In terms of these color coefficients, the complete amplitude for gg^ggA can be decom- 
posed into three separately gauge invariant sub-amplitudes 

^'' = E^^E'^f/- (2.17) 

1=1 f 

The sum over colors of the external gluons for the squared amplitude becomes 

J2\A^'\'=i2Af{AfyJ2c.c,, (2.18) 

color *>i=l color 
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where the color factors are given by 



Ci = CiCi = , (no sum. over i) , (2.19) 



color 



(iV2 _ 1) (3 - iV2) 



color 

Thus, one finally gets 

3 3 

J2\^''\'=CiJ2\-^f\' + ^'^ Yl -^fi-^TY- (2-21) 

color 1=1 i,j=l; ij^j 

3 Numerical implementation 

Analytic expressions for the amplitudes of the previous chapter were implemented in 
the Fortran program VBFNLO [IDl HZ]. The tensor reduction of the loop contributions 
up to boxes is performed via Passarino-Veltman reduction [11]. Additionally we avoid the 
explicit calculation of the inverse of the Gram matrix by solving system of linear equa- 
tions, which is numerically more stable close to the singular points. For pentagons, we use 
the Denner-Dittmaier algorithm [18] which avoids the inversion of small Gram determi- 
nants, in particular for planar configurations of the Higgs and the two final state partons. 
The program was numerically tested in several ways. Besides usual gauge-invariance and 
Lorentz-invariance tests of the amplitudes, the different topologies were checked sepa- 
rately. The contraction of a triangle-tensor T^^^^{qi, q2,rnf) with gluon momentum has 
to vanish due to total antisymmetry of the Levi-Civita symbol 

?r^WM2(9l>?2,m/) = gf T^iM2(?l>52,m/) = . (3.1) 

Contracting with external gluon momenta, the tensor expressions of boxes and pentagons 
reduce to differences of triangles and boxes, respectively. With the tensor integrals as 
defined in the Appendix, the Ward identities for the boxes read 

(li'Bf,,^^f,,{qi,q2,q3,mf) = T^^^^{qi2,q3,mf) - T^^^^{q2, qs^ruf) , (3.2) 

?2'^MiM2M3(9l'92,g3,"^/) =T^it^3{(ll,(l23,mf) - T^,^3 (gi2, ^3, "^/) , (3.3) 

{qi,q23,mf) , (3.4) 

where the abbreviation qij = qi + qj has been used. Similarly, for the pentagons one finds 
{quq2,q3,qA,mf) = B (gi2,g3,g4,m/) - B (g2,g3,g4,m/) , (3.5) 

q2^ B^l^J.2^i3^J.4,iq^y ?2) qs: q^y iTT-f) = -S^j^3^^ to,?4,"^/) - 5miM3M4(?12,93,?4,"2/) , (3.6) 
{qi,q2,q3,q4,'mf) = B^,^^^^^{qi,q2,q34,rnf) - B^^^,^^^{qi,q23,q4,mf) , (3.7) 
gr^'MiM2M3M4(9i,92,g3,g4,m/) = B^,^^^^.^{qi,q2,q3,mf) - B^^^^^^{qi,q2,q3A,mf) . (3.8) 
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These relationships were tested numerically and they, typically, are satisfied at the 10~^ 
level when using Denner-Dittmaier reduction for the tensor integrals. In addition, one can 
perform a QED-check for the pentagons. Replacing gluons by photons and considering 
the process '-ff^'j'jA, diagrams with three- and four-gluon- vert ices vanish, because these 
structures are not available in an Abelian theory. The amplitude is simply given by the 
sum of all pentagon graphs, without color factors. When contracting with an external 
gauge boson momentum one obtains zero, since boxes are not allowed for photons, by 
Furry's theorem. Our amplitudes pass this test as well. 

To check the full scattering amplitudes, one can make use of the heavy-top effective 
Lagrangian for a SM strength Yukawa coupling, 

Cfs = ^G^,,,,G^''''A with G-^mM^ = le/^iM.M3M4(ja^^^ _ (3.9) 

As rrit becomes large, the results calculated with full fermion loops must approach the 
approximate ones derived from the effective Lagrangian. This check was performed with 
rrit = 5000 GeV, and cross sections agreed well for Higgs boson masses in the range 100 
GeV< rriA < 1000 GeV. In production runs, we put a cut in the routines for the determi- 
nation of the tensor integral coefficients of the C and D functions such that the complete 
amplitude is set to zero when small Gram determinants appear. We have checked that 
the result and plots do not depend on this cut for a broad range of values. Finally, 
the amplitudes for all three subprocesses were recalculated using the FeynCalc/FormCalc 
framework [15]. For a selection of randomly generated phase space points, the two inde- 
pendent calculations yield agreement at least at the 10~^ level. 



4 Applications to LHC physics 

The numerical analysis of the $ + 2 jet cross section was performed with a parton 
level Monte Carlo program in the VBFNLO framework [T7], using the CTEQ6L1 [19] set 
for parton-distribution functions. In order to prevent soft or collinear divergencies in the 
cross sections, a minimal set of acceptance cuts has to be introduced. Following Ref. P, 
we impose 

PTj > 20 GeV , \r]j\ < 4.5 , % > 0.6 , (4.1) 

where ptj is the transverse momentum of a final state parton and Rjj describes the 
separation of the two partons in the pseudo-rapidity versus azimuthal-angle plane 

Rn = 4^?^, . (4-2) 

with Ar/jj = \riji — rij2\ and (pjj = 4'ji — (f)j2- These cuts anticipate LHC detector capabilities 
and jet finding algorithms and will be called "inclusive cuts" (IC) in the following. Unless 
specified otherwise, the factorization scale is set to 

Ai/ = y/PTl PT2 (4.3) 
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Figure 2: A+2 jet cross section of the individual contributions of the subprocesses 
quark-quark, quark-gluon and gluon-gluon scattering for tan (3 = 1 as a function of 
the pseudo-scalar Higgs boson mass, niA- Here, the inclusive cuts (IC) of Eq. (|4.1|) 
were applied. 



while the renormahzation scale is fixed by setting [T] 

= as{PTi)as{PT2)oil{mA) ■ (4.4) 

We use one-loop as running with as{Mz) = 0.13. All our results below contain the 
contributions from the full top- and bottom-quark loops. For the top-quark mass we use 
TUt = 173.1 GeV. In the case of bottom-loops, running Yukawa coupling and propagator 
mass are taken into account, with the Higgs-mass as reference scale. Within the Higgs- 
mass range of 100 to 600 GeV, the bottom-quark mass is 33 to 42 % smaller than the pole 
mass of 4.855 GeV. The evolution of rrib up to a reference scale /i can be expressed as 

mb(/i) = m5 (mb) '^^"^/^l^^^. , (4.5) 

with mi, {nih) = 4.2 GeV, as derived from the relation between pole mass and MS-bar 
mass. For the coefficient function c, the five fiavor approximation [201 EI] within the mass 
range < < m^, 

/23 \ 23 

c(x)=i—xj [l + 1.17549 x + 1.50071x2 + 0.172478 x^] , (4.6) 

is used. Further evolution of Mb to a renormahzation scale fi > mt can be performed 
safely within the five fiavor approximation, because the deviation to the six fiavor scheme 
is less than 1% for /i < 600 GeV. 

Contributions of individual subprocess categories to the total cross section for tan /3 = 
1 are shown in Fig. [2] as a function of the Higgs boson mass, niA- Here, the minimal cuts 
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Figure 3: A + 2 jet cross section as a function of the pseudo-scalar Higgs boson 
mass, iJiA for different values o/taii/3. The inclusive cuts (IC) of Eq. (j4.1|) are 
applied 

of Eq. fl4.ip were used. The cross sections for processes involving gluons (quark-gluon or 
gluon-gluon amplitudes) exceed the quark-quark scattering contributions by more than 
one order of magnitude. The dependence of the full cross section, with top- and 
bottom-quark interference, is given in the left panel of Fig. O for a selection of tan (3 
values. For small tan/3, amplitudes with a top-quark loop dominate over bottom-quark 
loop mediated contributions. The striking peak arises due to threshold enhancement near 
irtA ~ 2 mt, whereas for bottom-quark loop dominated processes the peak would appear 
well below the Higgs mass range shown. 

For low m^, the minimal cross section is obtained near tan P ^ 7, when ht ~ hf, (see 
Eq. f l2.3p ) and both Yukawa couplings are suppressed compared to hf'^^. For large tan/3, 
e.g. tan/3 = 50 in Fig. El^a), the bottom-quark loops dominate. However, they lead to 
a much more rapid decrease of the cross section with rising because the suppression 
scale of the loops is now set by the heavy Higgs boson mass instead of the quark mass. The 
reduced importance of the bottom-quark loops at large implies that equality of the 
top and bottom contributions and, thereby, the minimum of the production cross section 
is reached at increasingly larger tan /3 as ttia is increased. This effect is demonstrated in 
the right panel of Fig. [3l 

The 75-matrix in the Dirac-trace of the quark loops leads to a new tensor structure and 
to a normalization of the loops that, for equal Yukawa couplings, induces a (3/2)^ = 2.25 
times larger Ajj than Hjj cross section. This enhancement is shown in Fig. H] and 
is also apparent in the effective Lagrangian of Eq. fll.ip . where the coefficient of the 
CP-odd Agg coupling exceeds that of the Hgg coupling by a factor 3/2. This effective 
Lagrangian provides a good approximation to the total $jj cross sections up to Higgs- 
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Figure 4: Comparison of cross sections for the CV-odd and CV-even Higgs for 
the full loop calculation and within the effective theory. The inclusive cuts (IC) of 
Eq. (j4.1|) are applied. 



masses of ~ 160 GeV and for small transverse momenta, ptj ^ 'm-t- In this region, 
the effective Lagrangian approximation can be used as a numerically fast alternative for 
phenomenological studies [1]. 

The smaller quark mass in the bottom loops also has a pronounced effect on the 
transverse momentum distribution of the accompanying jets: for pTj ^ rrib the large scale 
of the kinematic invariants leads to an additional suppression of the bottom induced sub- 
amplitudes compared to the heavy quark effective theory. This effect is clearly visible in 
Figs. E] and El where the transverse-momentum distributions of the softer and the harder 
of the two jets are shown for pseudoscalar Higgs masses rriA = 120, 200 and 400 GeV for 
tan/3 = 1, 7, 30. For modest Higgs mass values, both distributions fall more steeply for 
large tan/3. At large values of rriA, the Higgs boson mass sets the scale for the fermion 
loops which, in the rriA = 400 GeV panels of Figs. Eland El leads to pt distributions which 
are approximately equal for the top- or bottom-quark dominated loops. 

The azimuthal angle between the more forward and the more backward of the two 
tagging jets, = (j)jF — (f>jB, provides a sensitive probe for the CP-character of the Higgs 
couplings to the quarks [221 EH EE]. As shown in the left panel of Fig. |H] for a heavy 
quark in the loop, the CP-even Hqq coupling produces a minimum for (pjj = ±90 degrees 
while a 75-induced CP-odd Aqq coupling leads to minima at (pjj = and ±180 degrees. 
The softening effects observed for the jet transverse momentum distribution then raise 
the question, to what extent the jet azimuthal angle correlations of the effective theory 
will get modified when bottom quark loops dominate. 

For the set of pseudoscalar Higgs masses and tan /3 values mentioned above, predictions 
for the normalized (/)jj-distributions are shown in Fig. [71 The calculation was carried out 
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Figure 5: Normalized transverse-momentum distributions of the softer jet in Ajj 
production at the LHC, for different tan j3 and Higgs-mass values. The inclusive 
selection cuts of Eq. ()4.ip are applied. 




Figure 6: Normalized transverse-momentum distributions of the harder jet in Ajj 
production at the LHC for different tan /3 and Higgs-mass values. The inclusive 
selection cuts of Eq. (14. ip are applied. 



with a modified set of cuts, however, which was shown in Ref. [25] to lead to a better 
sensitivity to the CP-structure of the Higgs couplings than the inclusive cuts. In Fig. [71 
we use 

PTj > 30 GeV , \r]j\ < 4.5 , % > 0.6 , Ar]jj = |r/ji - Vj2\ > ^ , (4.7) 

which we call ICphi cuts in the following. One finds that the characteristic structure of 
the (pjj distribution, dips at = and = ±180 degrees, remains for bottom quark 
dominated Ajj production, albeit at a quantitatively reduced level for > 2mq. For a 
relatively light pseudoscalar Higgs boson and large tan /3, the softer transverse momentum 
distribution of the Higgs leads to kinematical distortions of the distribution: at ^ 
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Figure 7: Distribution in the azimuthal- angle between the two final state jets for 
CV-odd Higgs boson production at different Higgs-masses and tan/3 values. The 
ICphi set of acceptance cuts (see Eq. jj-'^ ) is used for pp collisions at I4 TeV. 



the Higgs recoils against two jets and hence must have pth > 60 GeV, and this high pt- 
scale leads to an additional suppression as compared to the (pjj ~ ±180 degree case where 
transverse momentum balancing of the jets does allow pth = 0. 

For the effective theory of the large quark mass limit, it was observed that CP- violating 
effects due a mixture of CP-even and CP-odd couplings leads to a phase shift of the 
distribution compared to the CP-even case by an angle, a, which is given by the relative 
strength of the two couplings f2M [20] • Taking into account the relative enhancement by 
the factor 3/2 of the Agg coupling due to loop effects, the phase shift angle is given by 

tan« = -^ (4.8) 

2 Vq 

when heavy quark loops of a single flavor dominate. In order to test this effect for the 
case of a light quark, we show, in the right panel of Fig. |H1 the results for 

3_ „mi, 3_ rrit 

Vb = T,yb = - tan/3 — and yt = ^Vt = - cot (3 — , (4.9) 
2 V 2 V 

where a 45 degree phase shift is expected, with minima of the distribution at -45 and 
-1-135 degrees. This basic expectation is, indeed, confirmed by the detailed calculation. 
However, there are additional distortions of the azimuthal angle distributions which can 
again be explained by kinematical effects due to transverse momentum balancing of the 
two jets and the Higgs boson. 



5 Conclusions 

In this paper, we have presented the determination of quark mass effects on the cross 
section and on distributions for pseudoscalar Higgs production in association with two 
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Fi gure 8: (f>jj distributions for different Higgs sector scenarios: pure CP-odd 
or CV-even coupling in the effective Lagrangian limit (left panel) and the CV- 
violating case defined in Eq. 14-^1 (right panel). In the effective Lagrangian, the 
CV-odd coupling was matched to the CV-even coupling by a factor 2/3. For the 
CV -violating Higgs-sector, results are shown for different tan/3 values and with 
jet acceptance cuts tightened to A'qjj > 4 and pTj > 40 GeV as compared to the 
ICphi set K7\\ . 



final state partons. Our calculation for Ajj production complements the analogous one 
for a scalar Higgs, i.e. Hjj production as carried out in Ref. [1]. Qualitative features are 
quite similar for the two cases. Validity of the heavy quark mass approximation is found 
to be restricted to m$ < rUg and prj < rriq while large dijet invariant masses do not spoil 
the validity of the heavy quark limit. A pronounced difference between Ajj and Hjj 
production is observed in the azimuthal angle distribution between the two jets, which 
allows, in principle, to determine the CP-properties of the produced Higgs boson at the 
LHC [26]. 

Our analytical expressions have been implemented in the VBFNLO program [10] and 
are publicly available as a parton level Monte Carlo program. Even though the code 
must evaluate loop expressions up to pentagons, the calculation is leading order in the 
strong coupling constant since Higgs production in gluon fusion first appears at the one- 
loop level. As a leading order process, it has been provided with an interface in the 
Les Houches format [28] to run with parton shower programs, providing full particle, 
momentum and color flow information. The code allows to sum top- and bottom-quark 
induced contributions with arbitrary CP-violating couplings 

^Yukawa = QiVq + ''^'l^yq)q^ , (5.1) 

and, thus, is versatile enough for simulating the effects of general, CP-violating Higgs 
sectors at the LHC. 
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A Tensor structure of triangles 



qi 



q2 



^2 




qi 



q2 



^2 




Figure 9: Two three-point functions connected by charge conjugation. 

The generic three-point functions for triangle graphs with opposite loop momentum 
have the following expressions 



i fd^k 



4m f / ivr 



-tr 



— vp?^ (k + — ' [k + qyi)^ — 



1- 



m. 



(A.l) 



i r d^k 



Amt I ivr 



-tr 



X+rrif ^^^^ M+a^ + rrif ^^^^ _X±q^±rri^^5 
k"^ — {k + — m?j {k + 512)^ — Tn^f 



(A.2) 



where gi, g2 are outgoing momenta, qi2 = Qi + ^2 and the overall factor — i/4m/ cancels 
the explicit mass factor arising from the Dirac trace. Using the charge conjugation matrix 
C 

C'TmC'-' = -7J , C^,C-'=^J with (7 = 7%', C^ = l, (A.3) 
one can derive (Furry's theorem [12]) 

Tr'^'(gi,g2,m^) =Tf'^^(gi,g2,m^) =T^^''^(gi,g2,m^) • (A.4) 
Thus, the color structure simplifies to 

tr[t"it"2]Tf '^^(gi,g2,mg) +tr[t"2ri]Tf '^^(gi,g2,my) = ^2, m^) . (A.5) 

Evaluation of the Dirac trace yields 

T''^^^(gi,g2,m^) = £™^''^ Co{q,,q2,mf) . (A.6) 

Here, Cq denotes the scalar three-point function and e^^i^i'^ is the totally anti-symmetric 
tensor (Levi-Civita symbol), contracted with the gluon momenta qi and q2- 
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Figure 10: Two four-point functions connected by charge conjugation. 



B Tensor structure of boxes 



The analytic expressions for the charge-conjugated boxes are 



-75^1/^2/^3/ N 

Bi (gi,g2,g3,m/) 



Am 4 



-tr 



X+mf X+jfi + rrif 



'1 

m'j {k + — TTi'j 
/ (A; + gm)^ - mj 



X 7 



/«2 



(B.l) 



-^5/^1/^2/^3/ N 

^2 {Qi,Q2,q3,mf) 



-i /• d^A; 



4m f / ivr 



-tr 



k"^ — m?P {k + q-^Y — 'm-'j 



X 7 



M2. 



y^+^23 + "^/ ,,Mi ^+^123 + "•^Z 



7' 



(k + g23)^ - "^J (A; + ^123)^ - 



m. 



■7 



(B.2) 



where qi, q-2, and gs are outgoing momenta, = gj + g^ and gjjfc = g^ + gj + qu- From 
charge conjugation one gets 



^1 (gi,g2,g3,"^/) = (gi,g2,g3,m/) = S (gi, ga, gs, "^/) • (B.3) 

Further two permutations are obtained by cychc permutation of (1, 2, 3). The color struc- 
ture for the sum of the two diagrams is 

i,{t-H'^H'^^)B^,'^'^\q,,q2,q,,mf)+ii{e^^^^^ 



B 



■A»l/i2A'3 



(gi,g2,g3,my 



2^'^^'^^^ (91,92, g3,m^) 



(B.4) 
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The tensor structure of charge-conjugation related box diagrams, e.g. with gluon permu- 
tation (1, 2, 3), can be written as 



^MlM2M3^^^ _ I ^^43915293 gtJ-ltJ-2 _ ^/i2'Jl'?2'?3 gfJ-lfJ-'i _|_ ^^2/^39293 ^A*! _ c-M2A'3'?l'?3 ^Ml 



9i 



?2 



_j_ ^M2A'3l3ll?2 gMl _|_ ^Ml<?ig2g3 ^M2M3 _j_ ^MlM3<72g3 gM2 _ ^MlAt3gi<?2 gM2 _ ^^41^29293 gM3 



_j_ ^A'lM2gig3 gM3 _|_ ^fJ.ltJ.2fJ.3q3 gtJ.lfJ'2 _ ^A'lAt2A»3'?2 ^M1M3 _j_ ^MlM2A»3'?l ^M2A»3 _j_ ^MlA'S'?!*^ ^2 
+ gM2) ^ ^/.l/.2...2 [2 + ^^^3) + ^^M3] I ^^^^^^ 

- e^^'^^^'^^^ Co(gi + g2, ga, m/) - e'^^^'^'^^"?^ Co(gi, + gs, my) 
+ 26/^^/^3^2.3 D^^(g,,g2,g3,my) + 2£/'™'?3 D ^,{q^, q^, q^,m f) 

+ 2£'^™''^D^3(gi,g2,g3,my). (B.5) 



The Z)o and are four-point functions. Whereas the former denotes a scalar function, 
the latter can be expressed by the usual Passarino-Veltman decomposition [H] as 



/^^(gi, g2, gs, rn/) = gf I^n + g^ I^i2 + gs ^13 • 



(B.6) 



Note that after contraction with polarization vectors and quark current J21 , the 

expression ( IB.SP still contains terms with factors (ei ■ gi), (e2 ■ g2), (^21 ■ gs) even though 
they vanish, since gluon polarization vectors and momenta q^ are perpendicular to 
each other and the quark current J21 is conserved. However, these terms are important 
for numerical gauge checks, where the corresponding gluon polarization vector is replaced 
by its momentum. Since the virtual gluon has a non-zero gf, these terms give finite 
contributions. 



C Tensor structure of pentagons 
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Figure 11: Two five-point functions connected by charge conjugation. 



The two five-point functions connected by charge conjugation are defined by the ex- 
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pressions 



—1 

4m f 



-tr 



X- 



m 



k"^ — m 



2 ' 

/ 



{k + ^4)2 — rn?f 



X+XlA + 



{k + gi4)2 - m^r"*' 



{k + gi24)^ - m^.^ 



{k + gi234)^ 



ml 



r7 



(C.l) 



P2™^*(gi,92,g3,g4,m;) = ^ / ^tr 

^ X+^23 + "-^/ X+Yl23 + ^/ ,,^4 X+^1234 + "-^/ ,^,5 3) 

(fc + ^23)^ -m) ik + qi2zY -mj {k + 51234)^ - mj 

where qi, q2, qz and q^ are outgoing momenta (gjj = qi + qj and similarly for gj^fc and qijki)- 
The allowed color structures are given in pOj. The pentagon was reduced using the Denner- 
Dittmaier algorithm [18] and is available as a FORTRAN-subroutine in VBFNLO [IT]. 
The full analytic expression in terms of Passarino-Veltman reduction can be found in [T6] . 



A;^ — m^ (/c + ga)^ — m^ 
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